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ABSTRACT : 

This  report  is  a  companion  to  the  report,  "Analysis  of  Inventory 
Record  Accuracy"  (NPS55SoFs0071A)  and  supplements  that  report  by  pro- 
viding detailed  information  about  the  simulation  model  employed.   The 
model  is  described,  a  programming  guide  is  provided,  and  the  program 
and  its  flow  charts  are  given. 
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1.   INTRODUCTION 

This  report  is  a  companion  to  the  report,  "Analysis  of  Inventory 
Record  Accuracy"  ,  and  supplements  that  report  by  providing  detailed 
information  about  the  simulation  model  employed.   The  simulation  is 
useful  in  its  own  right. 

The  simulation  models  the  operations  of  a  single  manager,  single 
warehouse,  multi-item  inventory  system  under  continuous  review.   Special 
attention  is  given  to  the  introduction,  influence,  and  discovery  of  in- 
ventory record  errors.   The  simulation  can  be  useful  to  others  in  the 
study  of  the  effects  of  inventory  record  errors  or  the  efficiency  of 
various  physical  inventory  policies.   It  is  also  useful  in  the  evaluation 
of  transient  phenomena  associated  with  continuous  review  inventory  poli- 
cies; for  example,  the  quarterly  procurement  costs  for  a  given  set  of 
inventory  policies,  and  the  protection  or  number  of  stockouts  a  given 
item  actually  experiences.   In  any  event,  the  simulation  should  be  em- 
ployed to  evaluate  any  physical  inventory  policy  being  considered  for 
implementation  within  the  Navy  supply  system. 

The  next  section  describes  the  model  and  data  employed.   Chapter 
3   is  a  short  programmer's  guide  to  the  use  of  the  simulation  program. 
The  remaining  four  chapters  present  the  program  variable  labels,  the 
program  flow  charts,  the  program  itself,  and  the  item  data  employed  by 
the  authors. 
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2.   THE  MODEL  AND  DATA 

2.1   General  Approach 

The  system  modeled  is  that  of  a  single  manager,  single  warehouse, 
multi-item  inventory  system  corresponding  to  a  typical  Navy  Supply  Depot  ac- 
tivity.  The  simulation  includes  receipt  and  requisition  processing,  re- 
plenishment, and  physical  inventory  processes.   It  is  very  nearly  a  general 
purpose  inventory  system  model  with  a  multitude  of  uses.   The  exact  form 
of  the  model,  of  course,  reflects  a  great  deal  of  structure  concerned  with 
stock  record  error  generation  and  correction.   Every  attempt  has  been  made 
to  make  the  simulation  program,  written  in  FORTRAN  IV,  as  easy  to  under- 
stand as  possible.   The  labels  were  chosen  to  be  meaningful,  so  that  one 
may  view  the  program  as  the  manipulation  of  records  and  quantities  which 
are  real  entities  in  actual  practice. 

Two  inventory  options  are  available  within  the  simulation:   a 
complete  wall-to-wall,  with  the  interval  between  inventories  as  specified 
by  the  user;  and  a  selected  item  inventory  just  prior-to-buying,  which 
causes  a  scheduled*  inventory  to  be  conducted  on  an  item  just  before  a  buy 
is  made  on  that  item.   The  length  of  the  simulation  is  also  selected  by 
the  user.   Certain  statistics  are  computed  daily,  which  are  used  to  gener- 
ate annual  reports.   In  addition,  the  results  of  each  wall-to-wall  inventory 
are  printed,  in  conjunction  with  statistics  which  have  been  generated  during 
the  period  since  the  last  wall-to-wall  inventory.   Throughout  the  entire 


"The  following  terminology  is  employed  with  respect  to  selected 
item  physical  inventories.   A  spot  inventory  is  conducted  in  connection 
with  a  warehouse  refusal.   A  scheduled  inventory  is  performed  in  connection 
with  a  prior-to-buying  physical  inventory  policy. 


daily  routine,  as  well  as  during  the  inventories,  errors  are  being  generated 
based  on  data  presently  available.   At  the  end  of  the  simulation,  a  summary 
of  the  entire  run  is  printed,  which  can  be  used  to  evaluate  the  inventory 
option  employed. 

2.2  The  Data 

The  item  data  for  the  original  study  (1)  consists  of  information 
extracted  from  a  random  sample  of  505  items  taken  ten  months  after  a  wall- 
to-wall  inventory  was  conducted  at  NSC  Newport  in  1965.   Of  the  sample  of 
505  items,  all  187  DSA  items  were  chosen  for  the  present  simulation.   These 
187  items  represent  a  1%  random  sample  of  the  approximately  18,000  DSA 
items  stocked  at  NSC  Newport.   The  simulation  employs  a  stock  battery  of 
561  items  which  is  a  simple  triplication  of  the  basic  187  item  sample, 
though  the  number  of  items  in  the  stock  battery  is  arbitrary. 

The  item  characteristics  selected  for  input  to  the  simulation  are: 

-  serial  number,  arbitrary  (SERIAL)  ; 

-  unit  price  (PRICE) ; 

-  recorded  on-hand  quantity  (RECOH) ; 

-  actual  on-hand  quantity  (ACTOH) ; 

-  dues  outstanding  (DUES); 

-  reorder  quantity  (Q) ; 

-  reorder  point  (RP) ; 

-  mean  quarterly  demand  (DEAR) ; 

-  mean  procurement  lead  time  in  quarters  (LT) ;  and 

-  mean  absolute  deviation  of  quarterly  demand  (MAD). 


^refers  to  a  label  used  in  the  program  as  the  name  of  either  a 
datum  or  variable.   See  Chapter  4. 
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NAVSUPRANDFAC  (2)  wall-to-wall  physical  inventory  accuracy  data 
is  used  in  the  simulation  because  of  its  completeness.   The  data  indicate 
that  for  a  wall-to-wall  inventory,  the  count  is  correctly  made  92.9%  of 
the  time.   For  the  other  7.1%  of  the  time,  errors  are  assumed  to  be  normally 
distributed,  with  a  mean  error  quantity  of  zero,  and  a  standard  deviation 
derived  from  the  data  as  follows: 

ACTOH  Standard  Deviation 

1.87 

4.65 

2.20 

11.80 

When  an  error  is  generated,  the  magnitude  of  the  error  is  at  least  one  unit. 

Data  on  receipt  and  issue  errors  were  obtained  from  the  Quality 
Assurance  and  Internal  Review  Division  of  NSC  Oakland.   These  data  indicate 
that  95.99%  of  the  time,  the  quantity  ordered  equals  the  quantity  received, 
and  the  receipt  is  processed  correctly  (RECTOK) .   For  the  other  4.01%  of 
the  time,  errors  occur  as  follows: 

-  1.37%:  received  8%  more  than  ordered  (RECTEG) 

-  1.38%:  received  6%  less  than  ordered  (RECTEU) 

-  0.64%:  receipt  not  posted  (RECTNP) 

-  0.62%:  receipt  posted  to  wrong  stock  record  (RECTPW) . 

For  issue  processing,  the  issue  is  correctly  processed  97.74%  of  the  time 
(ISSOK).   For  the  other  2,26%  of  the  time,  errors  occur  as  follows: 


-  0.73%  overissue  (ISSEO) 

-  0.48%  overissue  by  7% 

-  0.18%  overissue  by  15% 

-  0.07%  overissue  by  30% 

-  0.75%  underissue  (ISSEU) 

-  0.50%  underissue  by  8% 

-  0.12%  underissue  by  20% 

-  0.13%  underissue  by  50% 

-  0.78%  issue  wrong  item  (ISSWID). 

Data  on  scheduled  inventories  were  also  provided  by  NSC  Oakland, 
and  indicated  that  the  scheduled  inventory  is  performed  correctly  96%  of 
the  time.   For  the  other  4%  of  the  time,  errors  occur  as  follows: 

-  1%:   errors  are  plus  or  minus  1 

-  1%:   errors  are  plus  or  minus  2 

-  2%:   errors  vary  from  5  to  100,  as  a  function  of  ACTOH. 

For  spot  inventories  the  assumption  is  made  that  an  accurate  re- 
conciliation is  made  97%  of  the  time  and  that  the  record  remains  unchanged 
3%  of  the  time.   This  is  a  crude  assumption,  but  fortunately  the  spot  in- 
ventory plays  a  very  minor  role  in  the  simulated  inventory  operations. 

Data  on  inventory  costs  were  provided  by  FMSO  (3) .   For  conducting 
a  wall-to-wall  inventory,  the  cost  is  estimated  to  be  $1.09  per  item.   For 
conducting  a  selected  item  inventory,  the  cost  is  estimated  to  be  $3.85. 
A  spot  inventory  was  estimated  to  cost  $3.92.   These  costs  include  recon- 
ciliation, in  addition  to  the  physical  counting  process. 


2.3   Underlying  Distributions 

DEAR  and  MAD  item  demand  data  from  FMSO  records  is  used  to  specify 
the  parameters  of  the  item  demand  generation  processes.   Realistic  simula- 
tion requires  that  requisitions  be  generated  over  time  in  some  fashion,  and 
that  the  total  number  of  items  requisitioned  per  quarter  approximate  DEAR. 
Eoth  the  time  between  successive  requisitions  on  an  item  and  the  size  of 
the  requisitions  must  be  specified.   To  assume  that  all  requisitions  are 
for  a  quantity  of  one  would  be  misleading  since  the  total  number  of  requisi- 
tions the  system  would  have  to  process  would  be  too  large.   Since  each 
requisition  and  the  attempted  issue  which  follows  from  the  requisition  can 
generate  errors,  it  is  important  that  a  variable  requisition  size  be  used. 

The  assumption  made  in  this  study  is  that  item  demand  follows  a 
"stuttering  Poisson"  distribution;  that  is,  that  the  time  between  successive 
demands  is  exponentially  distributed  and  that  the  requisition  size  has  a 
geometric  distribution  on  the  positive  integers.   The  method  of  moments  is 
used  to  estimate  the  parameters  of  the  stuttering  Poisson  from  DEAR  and  maD 
for  each  item.   Equations  for  these  estimates  are  developed  in  Appendix  B  of  (1) 

The  daily  operations  of  a  stock  point  are  dynamic,  and  any  model 
which  attempts  to  simulate  such  operations  must  provide  the  randomness  which 
is  needed.   In  the  above,  the  underlying  distributions  for  the  simulation 
were  stated,  with  no  indication  of  how  these  distributions  were  to  be  gener- 
ated.  In  order  to  provide  the  randomness  called  for,  and  to  meet  the  basic 
criteria  of  the  distributions  as  stated  above,  pseudo  random  numbers  are 
generated  which  are  inputs  to  subroutines  which  output  random  variables 
with  various  distributions  as  required.   The  stream  of  numbers  is  random 
in  that  it  meets  certain  statistical  tests  for  randomness,  and  is  pseudo 


in  that  any  given  stream  of  random  numbers  can  be  reproduced.   In  the  pres- 
ent simulation,  the  IBM  subprogram  RANDU  is  used  to  generate  three  differ- 
ent streams  of  random  numbers:   one  for  generating  demands,  one  for  gener- 
ating errors  of  various  types,  and  one  for  all  other  uses.   These  three 
different  streams  allow  comparisons  to   be  made  between  different  inventory 
policies.   It  is  assumed  that  a  valid  comparison  of  inventory  policies  is 
possible  only  when  the  stock  sample  being  used  faces  the  same  pattern  of 
denand  from  run  to  run.   This  scheme  also  provides  for  allowing  runs  to 
be  made  with  no  errors  at  all,  these  "clean"  runs  providing;  a  benchmark 
for  the  "dirty"  runs  (those  runs  in  which  errors  are  introduced) . 

2.4   The  Vector  Framework 

The  basic  framework  of  the  simulation  is  two  vectors:   a  stock 
record  or  item  vector,  and  a  buy  vector.    An  item  vector  has  various  com- 
ponents which  allow  the  item  to  be  identified,  and  which  allow  the  status 
of  the  item  to  be  maintained.   The  components  of  each  item  vector  are: 
SERIAL,  PRICE,  RECOH,  ACTOH,  DUES,  Q,  RP ,  DEAR,  LT,  and  MAD,  as  defined  in 
Section  2.2  and,  in  addition,  the  following: 

-  quantity  presently  backordered  (BO) ; 

-  total  dollar  value  of  buys  to  date  (CUMBUY) ; 

-  demand  parameters  (P,  NU) ; 

-  date  of  next  requisition  (NXTREQ) ; 

-  total  cumulative  demand  to  date  (CUMDMD) ;  and 

-  total  cumulative  backorder  days  to  date  for  this  year  (BODAYS) 


The  subscript   I   is  employed  throughout  the  simulation  to  refer 
to  the  Ith  item;  for  example,  if   I  =  25  ,  then  PRICE(I)  is  the  unit  price 
of  item  number  25. 

The  second  vector  is  the  buy  vector,  which  consists  of  the 
following  components: 

-  the  item  number  to  which  the  buy  applies  (INDEX); 

-  the  quantity  ordered  (ORDQN) ;  and 

-  the  due  date  of  the  material  (DUEDAT) . 

The  subscript   J   is  employed  throughout  the  simulation  to  refer  to  the 
Jth  buy;  for  example,  if   J  =  681  ,  then  INDEX(J)  holds  the  item  number 
for  which  the  681st  buy  was  made,  e.g.,  item  number  25.   The  use  of  these 
vectors  will  be  clarified  in  the  following  paragraphs. 

2.5   Program  Operations  and  Control 

The  reader's  attention  is  invited  to  the  system  flow  chart,  as 
well  as  the  detailed  flow  charts.   The  present  discussion  will  provide  a 
narrative  clarification  of  these  charts,  and  is  intended  to  provide  an 
appreciation  of  the  logical  construction  of  the  model. 

After  the  standard  initialization  procedures,  the  item  data  is 
read,  after  which  the  user,  assumed  to  be  controlling  the  program  through 
a  time-sharing  terminal,  is  queried  for  the  length  of  the  simulation,  the 
type  of  run  (clean  or  dirty),  the  demand  random  number  streams  initializer 
(INITRN),  and  the  inventory  option  desired.   The  value  of  INITRN  is  arbi- 
trary and  merely  provides  for  a  constant  demand  pattern  from  run  to  run 


if  desired.   The  simulation  is  now  independent  of  the  user,  and  proceeds 
as  follows.   Due  dates  are  computed  for  any  outstanding  dues  at  the  be- 
ginning of  the  simulation,  by  assuming  that  all  material  will  arrive  during 
the  first  sixty  days  of  the  run  in  accordance  with  a  uniform  distribution. 
Then  the  parameters  of  the  stuttering  Poisson  demand  generators  are  com- 
puted from  item  DEAR  and  MAD  data.   Dates  of  first  requisitions  for  each 
item  (NXTREQ)  are  computed  using  the  exponential  distribution,  and  all 
records  are  scanned  to  determine  if  buys  are  necessary.   (Conceptually, 
the  generation  of  dates  for  the  next  requisition  to  occur  is  a  type  of 
'event-store'  process,  in  which  the  occurrence  of  a  particular  type  of 
transaction  causes  the  generation  of  the  time  for  the  next  transaction  of 
the  same  type  to  occur.) 

The  decision  as  to  whether  to  make  a  buy  is  made  by  the  subroutine 
BUY  by  computing  the  inventory  position  (IP)  of  the  item,  defined  to  be  the 
recorded  on-hand  quantity  (RECOH)  plus  outstanding  dues  (DUES)  minus  back- 
orders  (BO) .   If  the  inventory  position  is  less  than  or  equal  to  the  reorder 
point  (RP),  an  order  is  generated  for  the  integer  multiple  of  the  order 
quantity  (0)  which  will  bring  the  inventory  position  up  to  a  point  between 
RP  and  RP  +  Q. 

Having  completed  the  initialization  procedures  for  day  one,  the  daily 
routine  begins.   Each  day  is  identified  by  an  integer  number;  the  present 
day  at  any  time  is  the  value  of  the  variable  TODAY.   Thus  the  simulation 
begins  with  TODAY  equal  to  one,  and  time-steps  through  to  TODAY  equal  to 
FINISH,  which  is  the  last  day  of  the  simulation. 

The  daily  routine  begins  with  receipt  processing,  which  consists 
of  scanning  the  list  of  outstanding  orders  to  see  if  any  dates  in  the  DUEDAT 


vector  match  the  date  in  TODAY.   If  there  is  a  match  for  a  particular  item, 
a  random  number  if  generated  to  determine  whether  the  receipt  will  be 
processed  correctly,  or  with  errors.   After  the  receipt  is  processed, 
a  check  is  made  to  determine  if  there  are  any  backorders  outstanding   for 
that  item.   If  so,  an  attempt  is  made  to  release  the  backorders.   (Note: 
even  though  requisitions  of  different  sizes  are  generated,  only  the  total 
quantity  of  each  requisition  backordered  is  recorded  in  BO.   No  provision 
is  made  to  distinguish  a  backorder  resulting  from  one  requisition  from  that 
resulting  from  any  other  requisition.   Thus,  the  simulation  can  only  keep 
statistics  on  the  number  of  requisitions  which  have  been  backordered.) 
Backorder  releases  are  very  similar  to  regular  issues,  in  which  the  record 
is  checked  to  determine  if  an  attempted  issue  should  be  made,  after  which 
the  actual  on-hand  quantity  is  checked  in  order  to  actually  effect  the 
issue.   In  the  case  of  a  backorder  release,  an  attempted  issue  quantity 
(ISSQN)  assumes  the  value  of  the  present  number  of  backorders  (BO),  and 
a  determination  if  made  of  whether  an  actual  issue  can  be  made.   If  so, 
the  subroutine  ISBOER  is  called  to  make  the  backorder  release  subject  to 
errors. 

After  each  issue  of  material  on  an  item,  whether  to  release  a 
backorder  or  to  satisfy  a  requisition,  a  check  of  the  new  inventory  position 
is  made.   If  the  new  inventory  position  is  at  or  below  the  reorder  point, 
a  buy  is  made  in  order-quantity  multiples  to  bring  the  inventory  position 
back  up  to  the  (RP,  RP  +  Q)  interval. 

After  the  receipt  and  backorder  processing  is  completed,  all  requi- 
sitions are  processed.   The  component  NXTREQ  is  checked  against  TODAY  to 
determine  if  there  is  to  be  a  requisition  today  for  the  particular  item  in 
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question.   If  so,  the  requisition  size  (REQSIZ)  is  generated  using  the 
geometric  distribution,  and  a  new  NXTREQ  is  generated  for  the  item  from 
its  exponential  inter-demand  time  distribution.   As  in  the  backorder  re- 
lease routine,  the  attempted  issue  quantity  (ISSQN)  is  determined  by  REQSIZ, 
and  both  the  recorded  on-hand  quantity  (RECOH)  and  the  actual  on-hand 
quantity  (ACTOH)  are  checked  to  see  if  the  issue  can  be  made.   If  necessary, 
a  warehouse  refusal  is  generated  which  results  in  a  spot  inventory  being 
taken,  subject  to  errors.   The  actual  issue  is  made  by  calling  the  subrou- 
tine for  making  issues  with  errors  (ISSERR),  the  flowchart  of  which  is  not 
included  since  it  is  so  similar  to  ISBOER. 

The  end  of  the  day  brings  the  daily  update  to  keep  track  of  such 
items  as  total  accumulated  unit  backorder  days  to  date  this  year  (BODAYS) , 
the  record  accuracy  at  the  end  of  the  day,  and  the  dollar  value  of  investment 
recorded  and  actually  held  on  hand.   If  today  is  the  day  for  a  wall-to-wall 
iiiventory,  the  subroutine  WALLOP  is  called.   If  today  is  the  end  of  a  year, 
the  annual  report  is  generated.   If  today  is  the  end  of  the  quarter,  the 
subroutine  QTR$  is  called  to  compute  the  dollar  value  of  demand  for  this 
quarter  and  the  dollar  value  of  buys  for  this  quarter.   If  today  is  the 
end  of  the  simulation,  the  summary  is  generated.   If  today  is  not  the  last 
day,  the  variable  TODAY  is  incremented,  and  the  daily  routine  begins  again. 

2.6   The  Output  '  :  ' 

The  original  purpose  of  the  simulation  was  to  determine  the 
effects  of  record  accuracy  upon  measures  of  effectiveness  and  costs  for 
the  system.   Accordingly,  the  output  of  the  simulation  was  designed  to 
allow  comparisons  over  time  of  certain  statistics,  as  well  as  to  allow 
evaluation  of  the  entire  run.   Table  1  is  a  typical  annual  report. 
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TABLE  1.    TYPICAL  ANNUAL  REPORT 


SUPPLY  PERFORMANCE  MEASURES; 

CUMREQ  =  2791 

COMPEL  =  2454  =  87.93  Per  Cent  of  CUMREQ 

PARTFL  =  316  =  11.32 

ACOMFL  =  2487  =  89.11 

APARFL  =  304  =  10.89 

ABOREL  =  123 

BOREL  =  123 

TOTAL  BODAYS  (EODTOT)  =   661227 

BUYS  =  722 

REFUSE  =  7 

ERROR  MEASURES: 

RECTOK  =  694  =  95.33  Per  Cent  of  Total  Receipts 

RECTEO  =  14  =   1.92 

RECTEU  =  10  =   1.37 

RECTNP  =  5  =   0.69 

RECTPW  =  5  =   0.69 

ISSOK  =  2850  =  98.04  Per  Cent  of  total  Issues  (Includes  BO  releases) 

ISSEO  =  22  =   0.76 

ISSEU  =  15  =   0.52 

ISSWID  =  20  =   0.65 

QTR  DEMAND $              BUY$ 

5  20710.69  20811.81 

6  20371.13  19452.00 

7  15642.81  15522.56 

8  1'774.37  17279.88 
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Each  measure  applies  only  to  the  year  immediately  preceding  the 
day  of  the  report.   The  measures  are  defined  as  follows: 

-  CUMREQ:   total  number  of  requisitions; 

-  COMPFL:   number  of  requisitions  completely  filled  on  demand; 

-  PARTFL:   number  of  requisitions  partially  filled  on  demand; 

-  ACOMFL:   number  of  attempted  complete  fills; 

-  APARFL:   number  of  attempted  partial  fills  (the  difference 
which  may  result  in  actual  versus  attempted  is  due  to  the 
effects  of  record  inaccuracies); 

-  ABOREL :   number  of  attempted  backorder  releases; 

-  BOREL:   number  of  actual  backorder  releases; 

-  TOTAL  BODAYS  (BODTOT) :   total  unit  backorder  days,  in  millions; 

-  BUYS:   number  of  buys;  and 

-  REFUSE:   number  of  warehouse  refusals. 

The  error  measures  are  the  same  as  those  stated  above  in  Section  2.2,  and 
represent  the  mean  accuracies  of  the  stochastic  processes  generating  re- 
ceipt and  issue  errors.   DEMAND$  and  BUY$  provide  the  dollar  values  of 
quarterly  demand  and  buys  for  the  quarter  indicated. 

When  the  wall-to-wall  inventory  option  is  selected,  a  report  of 
the  inventory  results  and  statistics  accumulated  since  the  last  wall-to- 
wall  inventory  are  printed.   Table   2   is  an  example  of  such  a  report. 
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TABLE   2.    TYPICAL  WALL-TO-WALL  INVENTORY  REPORT 

728  =  DAY  WALL-TO-WALL  INVENTORY  HELD 

89.48  =  PER  CENT  RECORDS  ACCURATE  JUST  PRIOR  TO  INVENTORY 
90.82  =  MEAN  PER  CENT  RECORD  ACCURACY  DURING  PERIOD 
88.77  =  MINIMUM  PER  CENT  RECORD  ACCURACY  DURING  PERIOD 
50456.00  =  MEAN  DOLLAR  VALUE  OF  RECOH  DURING  PERIOD 
45360.80  =  MINIMUM  DOLLAR  VALUE  OF  RECOH  DURING  PERIOD 
50372.13  =  MEAN  DOLLAR  VALUE  OF  ACTOH  DURING  PERIOD 
45386.03  =  MINIMUM  DOLLAR  VALUE  OF  ACTOH  DURING  PERIOD 

91.80  =  PER  CENT  RECORDS  ACCURATE  JUST  AFTER  INVENTORY 

Information  from  both  of  the  above  reports  is  accumulated  for 
the  summary  at  the  end  of  the  simulation.   The  summary  consists  of  two 
parts,  one  of  which  summarizes  information  collected  on  an  annual  basis, 
the  other  of  which  summarizes  information  related  to  the  wall-to-wall  in- 
ventory periods.   If  the  PBUY  inventory  option  is  selected,  the  periodic 
portion  of  the  report  contains  information  accumulated  annually.   Table 
3  is  a  typical  summary. 

The  annual  portion  of  the  summary  is  derived  from  the  individual 
annual  reports  which  are  generated  as  the  simulation  proceeds.   DIFF  is 
the  difference  between  COMPEL  and  PARTFL,  and  represents  those  requisitions 
which  resulted  in  a  backorder  for  the  full  amount  of  the  requisition  (a 
'no-fill').   The  periodic  portion  of  the  report  shows  the  inventory  period 
by  number,  the  mean  record  accuracy  for  the  period  (RECACC) ,  and  the  mean 
dollar  value  of  investment  during  the  period,  both  recorded  (RECOH$) 
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and  actual  (ACTOH$) .   This  latter  value  of  mean  actual  investment  repre- 
sents information  which  is  never  actually  available  to  the  stock  point 
manager,  and  which  will  be  used  to  cost  out  the  effects  of  errors  in  the 
system.   The  row  labelled  SUM  can  be  used  to  average  RECACC,  RECOH$,  and 
ACTOH$  for  any  number  of  periods  desired. 
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3.   PROGRAMMING  GUIDE 

In  addition  to  the  description  of  the  model,  it  is  deemed  necessary 
to  further  assist  the  potential  user  with  some  amplifying  remarks.   Some 
of  the  points  would  be  obvious  with  detailed  study  of  the  program,  other 
points  would  not  become  obvious  without  some  experience  in  running  the 
program. 

1.   The  item  input  data  consists  of  one  item  per  card,  formatted  as 
follows : 


CO 

CHARACTERISTIC 

FORMAT 

1-  3 

SERIAL 

13 

22-26 

PRICE 

F5.2 

28-31 

RECOH 

14 

33-36 

ACTOH 

14 

38-40 

DUES 

13 

42-45 

Q 

14 

47-50 

RP 

14 

52-57 

DEAR 

F6.2 

59-60 

LT 

F2.1 

62-66 

MAD 

F5.2 

The  serial  number  if  arbitrary;  there  is  enough  space  between  SERIAL  and 
PRICE  to  insert  a  stock  number,  if  desired. 

2.   The  simulation  is  currently  set  up  to  use  561  stock  records.   As 
such,  the  core  required  in  the  IBM  360/67  is  approximately  230K  bytes. 
This  storage  requirement  obviously  depends  upon  the  number  of  stock  rec- 
ords used.   In  addition,  it  is  quite  dependent  upon  the  activity  of  the 
items  in  that  the  buy  vectors  (INDEX,  ORDQN,  and  DUEDAT)  are  required  to 
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be  longer  for  a  more  active  item  sample,  since  the  number  of  buys  per  year 
will  be  greater  for  more  active  items  than  for  less  active  ones. 

3.  The  average  cpu  time  per  year  of  simulation  is  about  one  minute 
on  the  360/67,  if  3-4  years  or  more  are  run.   Actually,  the  time  required 
per  year  of  simulation  increases  from  year  to  year,  due  to  the  search 
procedure  employed  for  receipt  processing,  i.e.,  each  day,  the  entire 
DUEDAT  vector  is  searched  to  determine  the  receipts  due  in.   As  noted 
above,  the  DUEDAT  vector  becomes  longer  each  time  a  buy  is  generated. 

4.  The  simulation  files  are  set  up  as  follows: 

DATASET  NUMBER       FILE  NAME 

4  Item  input  data 

5  Terminal  input  data 

6  Terminal  output 
8  Printer  output 

Obviously,  these  files  can  be  set  up  to  suit  the  individual  user.   If 
the  program  is  to  be  run  as  a  batch  job,  the  input  from  the  terminal  can 
be  modified  for  card  input. 

5.  Following  is  an  example  of  the  terminal  I/O: 


EXECUTION  BEGINS. . . 

PRESENT  PARAMETERS  ARE: 

Si   INPUT 

INITRN=        100,RUN=  1,WALINT=   364.00000 

6.  END 

ANY  PARAMETER  CHANGES? 

n 

DURATION  IN  YEARS  OR  DAYS? 

y 

DURATION? 
2. 

INVENTORY  OPTION:   WALL  OR  PBUY? 

wall 

INVENTORY  INTERVAL  IN  DAYS? 

182. 
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The  values  of  INITRN,  RUN,  and  WALINT  are  preset,  and  may  be  left  un- 
changed if  desired.   A  value  of   1   for  RUN  produces  a  dirty  run:   a 
value  of   0  produces  a  clean  run.   WALINT  is  set  to  364,  in  case  the 
PBUY  option  is  selected,  so  that  annual  results  will  be  printed.   The 
duration  may  be  specified  in  either  years  or  days.   If  the  duration  is 
less  than  one  year,  an  abort  will  occur  when  the  program  reaches  the 
point  at  which  the  summary  statistics  are  to  be  printed,  since  the  annual 
portion  of  the  summary  statistics  does  not  include  the  first  year  in  the 
computation  of  the  means.   The  option  of  selecting  the  duration  in  days 
may  be  helpful  when  the  present  status  of  each  stock  record  at  the  end 
of  the  simulation  is  desired.   For  example,  one  might  wish  to  see 
specifically  what  occurred  on  day  722.   He  would  then  run  the  simulation 
for  721  days,  and  subsequently,  for  722  days,  electing  to  have  the  stock 
battery  vectors  printed  at  the  end  of  each  run. 

6.  E8  is  a  local  subroutine  which  allows  efficient  dumping  of  output 
periodically  onto  the  output  spool  under  CP/CMS,  which  is  the  present 
time-sharing  package  under  which  the  simulation  has  been  run.   The  purpose 
of  E8  is  to  avoid  filling  up  the  user  disk  file  with  output,  thus  causing 
an  abort  in  the  program.   The  user  should  be  aware  that  the  simulation 
can  produce  a  great  quantity  of  output  if  desired,  and  should  thus  ensure 
that  sufficient  precautions  are  taken  to  handle  the  output. 

7.  The  user  will  notice  that  in  many  arithmetic  statements,  .5  is 
added  to  the  right  hand  side.   This  produces  an  upward  rounding  off 
before  truncation,  since  many  quantities  dealt  with  in  the  simulation 
are  small  integer  numbers. 
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It  is  apparent  that  the  program  will  require  some  modifications  prior 
to  being  run  on  any  machine  other  than  that  for  which  it  was  designed. 
For  example,  the  instruction  to  call  E8  must  be  deleted.   The  dataset 
numbers  will  probably  have  to  be  changed,  etc.   The  authors  assume  that 
the  user  is  adept  at  programming,  however,  and  feel  that  the  changes 
which  will  be  required  are  nominal. 
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A.   SIMULATION  PROGRAM  LABELS 

ABOREL:  Number  of  attempted  backorder  releases  for  year. 

ACOMFL:  Number  of  attempted  completely  filled  requisitions  for  year. 

APARFL:  Number  of  attempted  partially  filled  requisitions  for  year. 

ACTOH(I):  Actual  on-hand  quantity  for  Ith  item. 

ACTOH$(I):  Dollar  value  of  actual  on-hand  inventory  for  the  Ith  day. 

B0(I):  Amount  of  material  backordered  for  Ith  item. 

BODAYS(I):  Number  of  backorder  days  for  Ith  item  for  year. 

BOREL:  Number  of  actual  backorder  releases  for  year. 

BODTOT:  Total  unit  backorder  days  at  end  of  year  for  that  year. 

BUY:         Subroutine  for  checking  inventory  position  on  an  item,  and 
initiating  a  buy  if  necessary. 

BUYS:  Number  of  buys  for  a  year. 

BUY$(I):  Dollar  value  of  buys  for  Ith  quarter. 

COMPEL:  Number  of  actual  completely  filled  requisitions  for  year. 

CUMREQ:  Number  of  requisitions  for  year. 

DBAR(I):  Mean  quarterly  demand  on  Ith  item  (data). 

DUEDAT(I):    Due  date  of  material  ordered  on  Ith  buy  for  the  item  number 
contained  in  INDEX (I ) . 

DUES(I):      Amount  of  material  due  in  for  the  Ith  item. 

FINISH:      Day  on  which  simulation  terminates. 

IFLAG:       Flag  generated  during  requisition  processing  which  determines 
for  ISSERR  subroutine  whether  ISSQN  should  be  COMPEL  or  a 
PARTFL . 

INDEX (I):     Contains  item  number  of  Ith  buy. 

INVOP:       Inventory  option  (WALL  of  PBUY) . 

ISBOER:      Subroutine  for  making  backorder  releases  with  errors. 
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ISSERR:  Subroutine  for  making  requisition  issue  with  errors. 

ISSOK:  Number  of  requisitions  issued  without  error  for  year. 

ISSEO:  Number  of  requisitions  overissued  for  year. 

ISSEU:  Number  of  requisitions  underissued  for  year. 

ISSWID:       Number  of  requisitions  in  which  wrong  stock  number  was 
issued  for  year. 

ISSQN:        Issue  quantity;  used  in  both  issue  and  backorder  release 
subroutines . 

LT(I):        Mean  procurement  lead  time  for  Ith  item  (data). 

MAD(I):       Mean  absolute  deviation  of  quarterly  demand  for  Ith  item  (data) 

NEWREC :       Item  number  of  new  record  chosen  at  random  by  subroutine 
NUREC. 

NU(I):        Parameter  for  exponential  demand  distribution  for  Ith  item 
(computed  from  data). 

NUREC:        Subroutine  for  randomly  selecting  record  in  same  price 
range  as  record  under  consideration. 

NXTREQ(I):     Day  of  next  requisition  on  Ith  item  (recomputed  every  time 
a  requisition  is  received  on  Ith  item). 

ORDQN(I):     Quantity  of  material  ordered  on  Ith  buy  for  item  number  con- 
tained in  INDEX (I) . 

P(I):         Parameter  of  geometric  distribution  employed  in  determination 
of  requisition  size  (REQSIZ)  on  Ith  item  (computed  from  data). 

PARTFL:  Number  of  actual  partially  filled  requisitions  for  year. 

PBUY:  Prior-to-buy  inventory  option. 

PRICE(I):  Unit  price  of  Ith  item  (data). 

RECACC(I):  Proportion  of  records  accurate  at  end  of  Ith  day. 

RECOH(I):  Recorded  on-hand  quantity  for  Ith  item. 

RECOH$(I):  Dollar  value  of  recorded  on-hand  inventory  for  the  Ith  day. 

RECTOK:  Number  of  receipts  processed  without  error  for  year. 

RECTEO:       Number  of  receipts  for  year  with  quantity  actually  received 
greater  than  quantity  ordered. 
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RECTEU:       Number  of  receipts  processed  for  year  with  quantity  actually 
received  less  than  quantity  ordered. 

RECTNP:       Number  of  receipts  for  year  with  no  posting  to  recorded  on- 
hand  and  dues  fields  of  records o 

RECTPW:       Number  of  receipts  for  year  with  quantity  posted  to  recorded 
on-hand  field  of  randomly  selected  record,  using  NUREC 
subroutine . 

REFUSL:       Number  of  warehouse  refusals  for  year. 

REQSIZ:       Requisition  quantity  generated  from  geometric  distribution 
upon  receipt  of  a  requisition. 

RP(I):        Reorder  point  for  Ith  item  (data). 

SERIAL (I):     Stock  (item)  number  (sequential  from  1,  2,  ...)  of  Ith  item. 

SPOT:         Subroutine  for  conducting  spot  inventory,  called  every  time 
a  warehouse  refusal  is  generated,  or  with  PBUY  physical 
inventory . 

TODAY:        Current  date  (integer  number,  begins  with  1). 

WALL:         Wall-to-wall  inventory  option. 

WALINT:       Interval  between  wall-to-wall  inventories;  provided  by  user. 

WALLOP:       Subroutine  for  conducting  wall-to-wall  inventories  and  pro- 
ducing certain  statistics  pertinant  to  the  period  just 
preceding  the  inventory. 
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5.   FLOW  CHARTS  FOR  SELECTED  PORTIONS 
OF  THE  SIMULATION  PROGRAM 
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SYSTEtl  FLOW  CHART 


Initialize  arrays  and 
variables . 

'  "\ 

Read  in  data. 

Query  user  for  length 
of  simulation,  inven- 
tory option,  type  of 
run  (clean  or  dirty), 
and  random  number 
streams  initializer 
(INITRN) . 


Generate  due  dates  for 
dues  outstanding  on 
first  day. 


Compute  parameters  P 
and  NU  for  geometric 
and  exponential  dis- 
tributions for  each 
item. 


Generate  time  of 
first  requisition  for 
each  item. 


Generate  first  day 
buys  based  on 
initial  inventory 
position. 
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SYSTEM  FLOW  CHART  (continued) 


Process  receipts  and 
release  backorders. 


Process  requisitions. 
Generate  time  of  next 
requisition.   Buy  if 
necessary . 


Perform  daily  update 
for  statistical 
information. 


Perform  wall-to-wall 
inventory  if  required. 

Produce  annual  report 
if  required. 

Is  today  the  last  day? 

■0 


STOP     N 
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PROCESS  RECEIPTS  AND  RELEASE  BACKORDERS 


0 


Check  for  receipt  to- 
day on  Ith  item. 


Generate  random  num- 
ber to  determine  man- 
ner in  which  receipt 
will  be  processed. 


Y 

Increment  RECTOK. 
Post  ORDQN  to 
RECOH  and  ACTOH. 

Process  correctly? 

(95.99%) 

Quantity  error  over? 


Quantity  error  under? 


Quantity  not  posted? 


N   (0.62%) 


Increment  RECTPW. 
Post  ORDQN  to  ACTOH. 


Generate  random  stock 
number  in  same  price 
range  as  original: 
CALL  NUREC 


Post  ORDQN  to  RECOH 
of  randomly  selected 
record . 


(1.37%) 


(1.38%) 


(0.6A%) 


Increment  RECTEO. 

Post  ORDQN  to  RECOH. 

Post  1.08*  ORDQN  to  ACTOH. 


Increment  RECTEU . 
Post  ORDQN  to  RECOH. 
Post  .94*ORDQN  to  ACTOH, 


Increment  RECTNP. 
Post  ORDQN  to  ACTOH, 


0 
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PROCESS  RECEIPTS  AND  RELEASE  BACKORDERS  (continued) 


RECOH=0? 


B0=0 


■0 


Increment  ABOREL 


RECOH  less  than  BO? 


Set  ISSQN=RECOH 
Subtract  RECOH  from 
BO.   Set  RECOH=0. 


Increment  REFUSE. 


Set  ISSQN^BO 

Set  B0=0 

Subtract  ISSQN  from 

RECOH . 


ISSQN  less  than  or 
equal  to  ACTOH? 


Make  issue  with 
errors:  CALL  ISBOER. 


ACTOH=0? 


Add  ISSQN  to  BO 


Add  (ISSQN-ACTOH)  to 
BO.   Set  ACTOH=0. 
Increment  BOREL. 


Increment  BOREL, 


BUY  if  necessary: 
CALL  BUY. 


Perform  spot  inven- 
tory: CALL  SPOT. 


0 


0 


Buy  if  necessary; 
CALL  BUY. 


All  records  checked? 


0 


€) 


© 
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PROCESS  REQUISITIONS 


Requisition  on 
ith  item? 


Generate  REQSIZ  from 
geometric  distribution. 


Increment  CUMREQ 


REQSIZ  less  than  or 
equal  to  RECOH? 


Set  ISSQN=RECOH.  ; 
Set  RECOH=0.  Add  i 
(REQSIZ-ISSQN)  to  | 
BO.  Increment  ! 
APARFL  ' 


\  ISSQN  less  than  or 
I  equal  to  ACTOH? 


Set  IFLAG=0 


Make  issue  with 
errors:  CALL 
ISSERR. 


Buy  if  necessary; 
CALL  BUY. 


V 


Add  ISSQN  to  BO, 


I  Take  spot  inven- 
;  tory:  CALL  SPOT. 


iBuy  if  necessary:- 
icALL  BUY.        i: 


■0 


Subtract  REQSIZ  from 
RECOH.   Set  ISSQN  = 
REQSIZ.   Increment 
ACOMFL . 


ISSQN  less  than  or 
equal  to  ACTOH? 


Set  IFLAG=1. 


Make  issue 
with  errors: 
CALL  ISSERR. 


Increment  REFUSL. 


Buy  if  neces- 
sary: CALL  BUY, 


ACTOH=0? 


Increment  PARTFL. 
Add  (ISSQN-ACTOH) 
to  BO.   Set  ACTOH=0 


Take  spot  inven- 
tory:  CALL  SPOT. 


ij  Buy  if  necessary; 
■  iCALL  BUY. 


Q 


V 
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PROCESS  REQUISITIONS  (continued) 


1 


?y 


Generate  time  of 
next  requisition: 
NXTREQ. 


d) 


All  records  checked? 


CONTINUE:  daily  update 
is  next . 


G 
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SUBROUTINE  ISBOER:   Issue  (release)  backorders  with  errors 


Generate  RN 


Issue  correctly' 


(97.74%) 


N 


Overissue? 


(0.73%) 


N 


Underissue? 


X 


:  (0.75%) 


(0.78%) 


Increment  ISSWID. 


■<r 


Generate  new 
item  in  same 
price  range: 
CALL  NUREC. 


ISSQN  less  than  or 
equal  to  ACTOH  of 
new  record? 


Subtract  ISSQN  from 
ACTOH  of  new  record. 


f 


Increment  ISSOK. 
Subtract  ISSQN  from 
ACTOH . 


Increment  ISSEO, 


Over  by  7%  ? 


(0.48%) 


Subtract  1.07* 
ISSQN  from  ACTOH. 


N 

y 

Subtract  1.15* 
ISSQN  from  ACTOH. 

n-tTar-     1-"   1  ^"Z      *> 

uver 

uy     J.  J  /o  ; 

(0.18%) 

N 


(0.07%) 


Subtract  1.30* 
ISSQN  from  ACTOH. 


Increment  ISSEU. 


\exit7 

v 


Y 

Subtract  .92* 
ISSQN  from  ACTOH. 

Under  by  8%  ? 

(0.50%) 

Y 

N 

Subtract  .80* 
ISSQN  from  ACTOH. 



Under  by  20%  ? 

(0.12%) 

N 

Subtract  .50* 
ISSQN  from  ACTOH. 



(0.13%) 

vEXI 
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6.   THE  FORTRAN  IV  SIMULATION  PROGRAM 
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7.   ITEM  DATA  EMPLOYED  IN  THE  ORIGINAL  STUDY 

The  item  characteristics  (Q,  MAD,  etc)  are  define^  in  Section  2.2 
of  this  report. 
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